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Abstract 

A two-dimensional body moves through a rarefied medium; the 
collisions of the medium particles with the body are absolutely elastic. 
The body performs both translational and slow rotational motion. It 
is required to select the body, from a given class of bodies, such that 
the average force of resistance of the medium to its motion is maximal. 

There are presented numerical and analytical results concerning 
this problem. In particular, the maximum resistance in the class of 
bodies contained in a convex body K is proved to be 1.5 times re- 
sistance of K. The maximum is attained on a sequence of bodies 
with very complicated boundary. The numerical study was made for 
somewhat more restricted classes of bodies. The obtained values of 
resistance are slightly lower, but the boundary of obtained bodies is 
much simpler, as compared to the analytical solutions. 
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1 Introduction 



Consider a homogeneous medium of point particles at rest in Euclidean space 
M. d , and a body moving forward through this medium. The medium is highly 
rarefied, so that mutual interaction of the particles is neglected. The interac- 
tion of the particles with the body is absolutely elastic. It is required to find 
a shape of the body that minimizes or maximizes resistance of the medium 
to its motion. 

When thinking of this kind of problems, one can have in mind an artificial 
satellite of the Earth moving on a relatively low (say, 100 -j- 200 km) orbit; 
one has to minimize resistance of the rest of atmosphere, or a solar sail; then 
one has to maximize the pressure of the flux of solar photons on the sail. 

In order to specify the problem, one has to describe the body's motion as 
well as to define the class of admissible bodies. In classes of convex bodies 
with translational motion, the minimization problem has been extensively 
studied. Newton pQ obtained the solution in the class of (three-dimensional) 
convex axially symmetric bodies of fixed length and width. Since 1993, 
there have been obtained many interesting results in classes of convex non- 
symmetric bodies [2]-[10j. 

Note that resistance can be written in the form R[f] = J J D (1+| V/| 2 ) -1 dx dy, 
where the function z = f(x,y), (x, y) G D describes the front part of a convex 
body. Here the (orthogonal) coordinates are chosen in such a way that the 
body's velocity equals (0, 0, 1). Thus, the problem amounts to minimization 
of the functional R[f]. 

This approach is not valid, as applied to nonconvex bodies. (A body is a 
bounded connected set with piecewise smooth boundary.) The reason is that 
the above formula is not true if particles can hit the body more than once. 
In fact, there is no simple analytic formula for resistance in the nonconvex 
case. 

In general it is not easy to calculate resistance even for (nonconvex) bodies 
with well-behaved boundary; however, usually one can construct a minimiz- 
ing body or a minimizing sequence of bodies. In the three-dimensional case 
infimum of resistance is typically equal to zero [TIj 112^1 . On the contrary, 
in the two-dimensional case infimum of resistance is positive and usually can 
be explicitly found [12]. In higher dimensions, d > 3, the answer is the same 

1 This result is obtained for classes of bodies of fixed length and width [12] and for 
classes of bodies containing a bounded set and being contained in its e-neighborhood [TT] . 
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as in the three-dimensional case: infimum of resistance equals zero. Note in 
passing that the problem of maximal resistance admits a very simple solution: 
the front part of the body's surface should be orthogonal to the direction of 
motion, or should be composed of pieces orthogonal to this direction. 

It is also interesting to consider rotational motion of the body. Imagine 
an artificial satellite without orientation control; one can expect that in the 
course of motion, it will perform a (perhaps very slow) rotation. The problem 
of minimal mean resistance for non-convex rotating bodies in two dimensions 
was considered in [13]; it was proved that the gain in resistance, as compared 
with the convex case, is smaller than 1.22%. 

Here we study the problem of maximal mean resistance for rotating bodies 
in two dimensions, d — 2. This problem is far from being trivial, contrary 
to the case of purely translational motion. To see it, consider a unit disk 
on the plane, which is moving forward and at the same time slowly (and 
uniformly) rotating. Denote the disk by K\. "Cut off" a small portion 
of the disk contained in the ^-neighborhood of dK\ (e <C 1); the resulting 
set B is such that B C K\ C N £ (B) (here Af e designates e-neighborhood). 

mi j_ • • i 1 ,1 r • , Resistancc(_B) 1 o 

lhe question is: how large can the increase of resistance Resis tance(Ki) " e - 
Some estimates can be made immediately. Firstly, it cannot exceed 1.5. 
This hypothetical maximal increase is achieved if the velocity of a reflected 
particle is always opposite to the incidence velocity, v + = — 1>; in this case 
the momentum transmitted by the particle to the body is maximal. Next, 
if the circumference dK\ is partitioned into several small arcs and each arc 
is substituted with a pair of legs of a right isosceles triangle contained in 
K\ (the resulting body is shown on Fig. (a) below) then resistance increases 
approximately a/2 times. More precisely, if the length of each arc is 2e 
then the resistance increase is y/2sme/e ~ \/2(l — e 2 /Q); for the proof see 
Appendix 1. 

Another example is the body obtained by making deep and narrow rect- 
angular "hollows" on the boundary of K\\ see Fig. (b). Put the depth of any 
hollows to be e, width e 2 , and the distance between neighboring hollows, e 3 . 
Then approximately one half of the particles incident on a hollow get out 
with the velocity opposite to the initial one, and the rest of the particles get 
out with the velocity symmetric to the initial one with respect to the smaller 
side of the rectangle. Resistance of the body is, approximately, the arith- 
metic mean of the disk resistance and the (hypothetical) maximal resistance 
(1.5 times the disk resistance); that is, the resistance increase is 1.25 + o(l), 



3 



e -> 0. 




(a) (b) 

In this paper, the resistance maximization problem is studied 

(A) in the class of planar sets of the form r < 1 — ef(6/e), < e C 1, 
/ > in polar coordinates r, 9 and 

(B) in the class of sets contained in a fixed two-dimensional convex body. 
These problems still have to be rigorously stated; it is made in Section 

2, and the corresponding reformulated problems take the names (A) and 
(B). The restricted problem (A) seemed to be more amenable to numerical 
study and was examined first. It was not completely solved; nevertheless, 
we present here some numerical results. We believe that they are of interest, 
since they allow one to reach values up to 1.446 (which is rather close to the 
upper bound of resistance) by using relatively simple geometric shapes. 

When working on the restricted problem, there was gained experience 
that eventually allowed one to solve the problem (B). The answer here is 1.5; 
the corresponding maximizing sets have much more complicated boundary 
than those used in (A). 

The paper is organized as follows. The mathematical formulation of the 
problems is given in Section 2. The problem (A) is studied numerically in 
Section 3, and the problem (B) is solved analytically in Section 4. The ob- 
tained results are discussed in Conclusions. Finally, some auxiliary formulas 
related to resistance of zigzag shapes are derived in Appendices 1 and 2. 

2 Statement of the problem 

Let B C M 2 be a bounded connected set with piecewise smooth boundary; 
consider the billiard in M?\B. Consider a billiard particle that initially moves 
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freely, then makes several (at least one) reflections at regular points of dB, 
and finally, moves freely again. Denote by convS the convex hull of B. 

The trajectory of the particle intersects <9(convI?) twice: when getting 
in the set convB and when getting out of it. Note that if the point of 
first intersection belongs to Iq := <9(conv.B) fl dB, then the two points of 
intersection coincide. Introduce the natural parametrization of <9(conv.B) 
by the parameter £ G [0, L], where L = |<9(convI?)| is length of the curve 
<9(convB). 

Let £ and £ + be the first and second intersection points, and let v and 
v + be the particle velocity at these points, respectively. Denote by (• , •) the 
scalar product. Let be the outward unit normal vector to <9(convI?) at 
the point corresponding to £. For a vector w such that (w,n^) > 0, let us 
agree that the angle between and w is counted from to w clockwise 
or counterclockwise; in the first case it is positive, and in the second case, 
negative. Thus, the angle varies in the interval [— 7r/2, 7r/2]. Denote by (p 
and ip + the angles the vectors and form with —v and v + , respectively. 
(Note that one always has = n^+.) Thereby, the one-to-one mapping 
T B : (<p, £) — > (</? + ,£ + ) is determined. It is defined and takes values on a full 
measure subset of [— 7r/2, tt/2] x [0, L]. Moreover, the following holds true. 

(i) Tb preserves the measure \i given by d[i(<p,£) = cos <p dipd^; 

(ii) T B l = T B . 

These relations follow from the measure preserving property and from time- 
reversibility of billiard dynamics; see [13] for more details. 
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Suppose now that the center of mass of the body moves forward at the 
velocity e 2 = (0, 1) and the body rotates around the center of mass with a 
small angular velocity u < 1. Thus, each individual particle interacts with 
the body in the same way as if it there were no rotation. Resistance of the 
medium is a periodic vector- valued function of time, Rs(t), with the period 
T — 2ir/u. In order to derive the formula for the mean value of resistance 
7Z(B) — f L -Rs(t) dt, consider a reference system moving forward at the 
velocity e 2 . In this reference system, the body rotates around a fixed point, 
and there is a flux of particles of velocity — e 2 incident on the body. Each 
particle transmits to the body a momentum proportional to v — v + , where 
v = — e 2 = (0, —1) is the initial velocity of the particle, v + = (sin(<^ + — 
</?), cos((/? + — <£>)) is its final velocity, tp = (poi+out, t is the moment of the first 
intersection of the particle with <9(conv.B), ip + = (pt B \(ip,£)] here £ means 
the point of the first intersection of the particle with <9(convi?) and 
is the first component of the mapping T B . The mean resistance 7Z{B) is the 
sum of all momenta transmitted to the body in a time interval of length T, 
divided by T, that is, 



K{B) = -cJ J (sin((pf B) (<p,€)-(p), l+cos(^ B) (^,0-^))^(^>,0; 

(1) 

the ratio c is proportional to the medium density. Later on we shall specify 
the appropriate value of c simplifying the subsequent formulas. 

Further, changing the variables (jp,£) i-> (jp,£) = (jP( B) (jp,€), -^(v^O) 
in the integral ([1]) and taking into account that the Jacobian related to this 
change is 1 (property (i)) and that <pt B \{<p,£) = V 9 (property (ii)), one gets 

TZ(B) = -cJ J (sm(ip- <pf B) (<p,£)), 1 +cos(ip- tf B) (<p,i)f) dn(tp,i). 

^ ° (2) 
Comparing and d2J) and taking into account that sine is odd, one concludes 
that the first component of 7Z(B) is zero, that is, 

n(B) = -cj (l + cos(^+ ) (^e)-^))d//(^0-e 2 . (3) 

Let us now reduce the formula ([3]) to the form more convenient for com- 
putation. The curve d(comrB) is the union of a finite or countable family of 
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sets J , h, I 2 , • • ., 

<9(convB) = Uj/j. 

Here I = d(convB) fl dB is the "convex part" of the boundary OB, and 
<9(convI?) \ dB is the union of open intervals Ji, / 2 , • • •■ Respectively, convi? 
is the union of a finite or countable family of sets Qq, &>i, ^2, • • •, 

convi? = Ujfij, 

where Qq = B and the sets Q\, f2 2 , ■ ■ ■ are connected components of convB\B 
( "cavities" on B) . The enumeration is chosen in such a way that Ii C dQi (see 
the figure above). Note that for some sets there may happen <9fi, C dB. 
These sets ("interior cavities" of B) have no influence on resistance and will 
be ignored in the sequel. 

Each interval Ij (i 7^ 0) corresponds to an interval Li (modulo L) on the 
parameter set [0, L], and Jo corresponds to [0, L] \ Uj^oA = : Lq. Let U be 
the length of Lj. The sets [—tt/2, it/2] x Li are invariant with respect to Tg; 
denote by Tj the restrictions of T B on these sets, and by <Pi((p, £), £i((f, £), 
the components of Tj. One easily sees that T is given by T Q (ip, ^) = (—</>, 0- 
The mean resistance 7Z(B) ([3]) is the sum of partial resistances TZi related to 
the ith cavity: 7Z(B) = ^T^i, where 

Hi = -c f [ (1 + cosfo+(p, - ■ e 2 . (4) 

J-n/2 J d 

In particular, 

/•tt/2 /• g 

1Z = — c / / (1 + cos 2ip) cos if dtp d£ ■ e 2 = — c/oe 2 . 
J -tt/2 J c 3 

Put c = 3/8, then one gets a slightly simplified relation: TZo = —lo e 2 . 

We say that a bounded set Q C M 2 with piecewise smooth boundary is a 
standard cavity if it contains the interval X := [0, 1] x {0} and is contained 
in the upper half-plane {(xi,x 2 ) : x 2 > 0}, that is, 

ICOC {{xi,x 2 ) : x 2 > 0}. 

Consider billiard in Q; suppose that a billiard particle starts from a point 
of X, and after several reflections from dQ \ X returns to X. Let (£, 0) and 
(sin 99, cos ip) be the initial location and initial velocity of the particle, and 
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denote by (£n(<£, an d ~( s i ia( Ptt{ l P,£),cos(Pn( ( P,0)i the fi na l point and 
final velocity. The so denned map (<£>,£) i— > (<£>n(y, £)> £^(^0) preserves 
the measure /i and is defined and takes values on a full measure subset of 
[-i/2, tt/2] x (0, 1], 
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1 

Designate 

J7 ( n ) = q / / C 1 + cos(^n(v,0 -^)) cos <pd£d<p; (5) 

The integrand in (jSJ) does not exceed 2, therefore JF(ft) < (3/8) Jq 1 2 cos tp d^dip = 
1.5. On the other hand, denoting ft e := 7r e ft, where ir e : (21,22) > (xi,ex 2 ), 
one has that fnAf^O ~¥ as £ ~* + > hence lim e ^ + .F(ft e ) = (3/8) f^^ 2 J 1 (l+ 
cos2ip) cos(pd!;dip = 1. It follows that 1 < sup n .F(ft) < 1.5. 

Each pair (ftj,ij), i = 1, 2, . . . can be reduced, by a similarity transfor- 
mation and a translation, to the form (ftj,X), where ft, is a standard cavity. 
Denote Aj = /j/|<9(conv£>)|; one has ^ Aj = 1 and IZi = — AiJF(ftj), therefore 

= -|<9(conv5)| ■ (A + ^ A^(ft t )) • e 2 . (6) 

Let the set B £ be given by r < l—ef(9/e) in polar coordinates r, 0, where 
/ is a 1-periodic continuous piecewise differentiable nonnegative function and 
s divides 2ir. Suppose that /(0) = f(l) = 0. As e — > 0, 1Z(B £ ) tends to 
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—2tt ■ !F(Qf), where 

Q f = {(x 1 ,x 2 ) : < x x < 1, < x 2 < f{x 1 )}. 

Thus, the problem (A) (which was not yet rigorously formulated) can be 
stated as follows: 

(A) Find sup^ FiQf) over all continuous piecewise differentiable nonneg- 
ative functions f : [0, 1] — > M + such that f(0) = f(l) = 0. 

The problem (B) reads as: find sup BcA - \TZ(B)\, where K C M. 2 is a 
convex bounded set with nonempty interior. In view of (J6j), it amounts to 
the problem 

(B) Find sup n jF(fi) over all standard cavities Q. 

Indeed, let Q n be a sequence of sets solving the problem (B); then a 
sequence of bodies B„ C K approximating K solves the problem (B), if all 
the cavities of B n are similar to Q n and the length of the convex part of B n 
tends to zero. Thus, one has sup BcE - \1Z(B)\ = \dK\ ■ sup^ ^(fi). 

3 Problem (A): numerical results 

Introduce the shorthand notation T[f\ := ^(Qf). Note that the functional 

T is continuous in the C topology: if f n — > f then J-[fn\ — > J~[f\] thus 
there exists a sequence of piecewise linear functions maximizing T . Therefore 
it seems natural to look for the maximum in classes of continuous functions / 
with piecewise constant derivative /'. We also examined classes of continuous 
functions with piecewise constant second derivative /". In the first case the 
graph of / is a broken line, and in the second, a curve composed of arcs of 
parabolas. 

In each numerical experiment there were made N1N2 trials with a billiard 
particle in Qf. Usually N\ and A2 were taken equal and varied from several 
hundreds to several thousands (up to 5000). Initially, the particle is located 
at 0) G J, where & = (i — l/2)/iVi, and has the velocity (smipj, cos^j), 
where ip$ = vr(j - ^-)/N 2 , i = 1, . . . , N u j = 1, . . . , iV 2 . Then the least 
time instant, when the particle gets into X again, is fixed, and the particle 
velocity = — (sin iff,, cosipf) just before this instant is registered. The 
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sum 

3 7T Nl N * 

' F= 8 WW ^ cos ^ i 1 + cos (^ ~ vtj)) 
1 2 »=i i=i 

is considered to be an approximation for the integral T[f}. 

The algorithm simulating the billiard dynamics, as well as numerical in- 
tegration, were implemented in Programming Language C. The precision 
accuracy achieved 10~ 6 ; it was controlled, firstly, by the differences between 
the successive approximations of T as N\ and N 2 increase, and secondly, by 
comparison with the analytic results. To maximize resistance, there were 
used optimization algorithms from the Genetic Algorithm and Direct Search 
Toolbox of the computational system MATLAB, version 7.2; these methods 
do not require any information of the derivatives of the objective function. 

The obtained results follows. 

1. In the class of two-segment broken lines 

f m _ J «£, if 0<£<£o 

I 0(1-0. if £o<£<i, 

where a > 0, (3 > 0, < £o < 1, «£o = 0(1 — £o), the maximum of T 
equals 1.42621 and is achieved at a = (3 = «o ~ 1.12. Then £ = 0.5; the 
corresponding set £lf , is an isosceles triangle with the angle 83.6° at the 
top vertex. It is shown on Fig. (a), with the lateral sides drawn dashed. 
The function J-[f a>a ] oscillates and goes to 1.25 as«-> +oo (see Fig. (b)). 



1.4 
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2. In the class of two-segment piecewise quadratic functions 



/aia 2 ^i/3 2 (0 - 



where < £ < 1) a iQ + Piio = 0:2(1 — £o) 2 + $2(1 ~~ £0), the maximum of T 
is achieved at £ — 0-5, cti = a 2 = a ~ —0.486, Pi = j3 2 = ~ 1-361, and 
is equal to 1.43816. The corresponding set Qf aaf}(} is a curvilinear isosceles 
triangle; it is shown on Fig. (a) above with lateral sides drawn with solid lines. 
Its height is equal to the height of the optimal triangle from the item 1; so 
to say, this triangle is obtained from the previous one by a slight "bending 
outwards" its lateral sides. 

3. In the class of broken lines with many segments the simulations be- 
come more cumbersome. Let x° = (x^x®) = (0,0), x 1 = (x\,x\), . . . ,x m = 
(x™,x™) = (1,0) be the vertices of the broken line, with m being the num- 
ber of segments. Experiments with relatively small m (m < 5) showed that 
making the broken line symmetric with respect to the vertical line x\ = 1/2 
and taking the values x\, i — 1, . . . , m — 1 equally spaced in [0, 1] favor some 
moderate increase of resistance. Therefore the posterior study was restricted 
to symmetric broken lines with x\ = i/m, and thus the number of effective 
parameters was reduced almost four times: from 2(m — 1) to [ttJ- This deci- 
sion allowed one to take a relatively large number of segments, 1 < m < 18. 

There were found many "zigzag" shapes providing nearly maximal values 
T ~ 1.446...; some typical shapes are shown on the figures (a,b,c) below. 
The greatest found value of T corresponds to the shape shown on Fig. (b) 
(m = 10) and equals 1.446227. 




(a) m — 6 




(b) m = 10 



(c) m — 18 
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In order to verify the simulation results, as well as to find out the com- 
mon analytic form of the obtained curves, we examined the following broken 
line. Fix \1/ e [0, 7r/2] and consider the arc of angular size 2\I/ contained 
in the upper half-plane %2 > 0, with the endpoints (0,0) and (1,0). Let m 
be even. Mark the points x° = {x^x®) = (0,0), x 2 = (xl,x?,), . . . , x 2t = 
(xf,xf), ...,x m = {xT,x%) = (1,0) on the arc, with = x\ < x\ < 
. . . < x™ = 1. Let 5 be the maximum of values xf — x 2 * -2 . We say 
that a right triangle ABC is canonical, if the vertex B is situated above 
the hypotenuse AC and the median drawn from B to AC is vertical. For 
i = l,...,m/2, draw the canonical triangle Ai 2l_2 i 2i_1 i 21 with the hy- 
potenuse [x 2i ~ 2 , x 21 }. Thus, one has xf^ 1 = \ (xf' 2 + xf), x 2 ^ 1 = \ (x 2 2 l - 2 + 

xf) + \y/{xf - xf~ 2 ) 2 + (xf - xf' 2 ) 2 . The broken line x ^ 1 . . .x m ~ 1 x m 
composed of legs of all triangles obtained this way will also be called canon- 
ical (the graphs shown on the figures (a,b,c) are good approximations for 
canonical lines with m — 6, 10, 18). If ^ is fixed and S goes to zero, the 
corresponding value of T tends to 

1 , 2v/2 sin | - 2 sin 4 | - * 

T m = 1 + - sin 2 * + — 2 - 2 - 7 

6 sin^ 

(here ^ is expressed in radians). The proof of this convergence is put in 
Appendix 2. The maximal value of is achieved at \l/o ~ 0.6835 ~ 39.16° 
and is equal to ^(^ ) = 1.445209. 

Let F m {^) be the values of T related to the canonical lines with x\ = i/m, 
i = 0, . . . , m. These values were numerically calculated for various values of 
\1/ and for m = 6, 10, 18. The resulting functions jF m (^) and the function 
JF(^) (solid line) are shown on the figure below. 
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Note in passing that in the limit \P — ► the triangles of the corresponding 
canonical line approach a right isosceles triangle, therefore the corresponding 
value of T tends to \/2; that is, .F(O) = v2. 

4. In the class of piecewise quadratic functions, with < 18 segments, 
the greatest found value of T is 1.44772; this is the maximal value found 
numerically up to the moment. The corresponding curve is a 14-segment 
"zigzag" line with slightly concave segments; its visual perception is the 
same as of the curves shown above on Figs.(a,b,c). 




Thus, one has the following estimation for the Problem (A): 

1.44772 < supJ r (fi / ) < 1.5. 

/ 

5. We examined numerically shapes formed by infinitely small canon- 
ical triangles. In these experiments we had to substitute the usual billiard 
dynamics with the pseudo-billiard one described in Appendix 2. Also, we 
tried shapes formed by various kinds of non-canonical triangles. However, 
we could not increase resistance this way. 
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4 Solution of Problem (B) 

Consider the standard cavity Q £ = Q! e U ft", where Q' £ = [0, 1] x [0, e] is 
a rectangle and Q" is the upper semi-ellipse with the foci F\ = (0,e) and 
F 2 = (l,e) and with the major semi-axis of length 1/e. This figure is a 
mushroom, with the stem and the cap Q". Note that mushroom was first 
proposed by Bunimovich as an example of billiard with divided phase space 

H3 




Consider the billiard particles starting from I = [0, 1] x {0}. These 
particles, except for a small part of them, make a unique reflection from the 
elliptical arc and then return to X, the angle between the initial and final 
velocity being less that 2 arctan(e/2). The rest of the particles, i.e. those that 
make reflections from the vertical sides of the stem, have the total measure 
0(e). Indeed, one easily calculates that the measure of the particles having 
the first reflection from the stem is 

2- d£ cos (p dcp = 2(1 + e- Vl+e 2 ). 

JO </arctan(£/e) 

The measure of the particles having at least one reflection from a stem size 
is at most double this value. 

Taking into account that the total measure of all particles incident on I 
is Jo cos ^ d(,d(p = 2 and using (jSJ) , one gets 

F(^e) > o ( 2 _ 4 ( X + £ ~ Vl +£ 2 )) f 1 + cos(2 arctan J)) = 1.5 + 0(e). 
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Thus, problem (B) is solved: sup^jFf^) = 1.5. 

Now let K be a convex bounded body with nonempty interior. Approx- 
imate it by a convex polygon K £ C K, e > such that \dK\ — \dK e \ < e. 
To each side a, b, ... of K £ assign a rectangle R a , Rb, . . . such that one side 
of the rectangle (denote it by a', b',..., respectively) belongs to a, 
etc; all the rectangles belong to K £ and do not mutually intersect; and the 
common length of the part of perimeter of K £ not occupied by a', b', . . ., is 
less than e. 




On each rectangle plant out a "seedlings of mushrooms", as shown on 
the figure (a) below for R a . The sides of R a that do not belong to a are 
shown dashed. On Fig. (a), there is shown R a \ (union of mushrooms) , all the 
mushrooms being of equal size and similar to fl £ . The total length of the 
lower horizontal part (l.h.p.) of the boundary of the obtained figure is 1 — e 
times the length of the corresponding size of R a , that is, (1 — e)\a'\. 

Now, plant out a "seedlings" of "mushrooms of the second order" (see 
Fig. (b) below). Here the length of l.h.p. of the boundary is (1 — e) 2 |a'|. 



Ra 




(a) 
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Exxxa [Z^S^l Exxx^ Exl 

(b) 



Continuing this process, one finally obtains the figure R a such that the 
length of the l.h.p. of its boundary is less than e\a'\. Similarly, one obtains 
the figures R b , . . .. 

Put B £ = (K £ \ (R a U R b U . . .)) U (Ra U R b U . . .)■ All the cavities of 
5 e and similar to Q e , the length of the convex part of dB £ tends to zero 
and |<9(conv.B e )| — > \dK\ as e — > 0. (Recall that the convex part of dB £ is 
dB £ nd(convB £ )). Therefore lim £ ^ \TZ(B £ ) \ = 1.5\K(K)\; this solves the 
maximal resistance problem (B). 

5 Conclusions 

There are many technical devices utilizing wind pressure force: for example, 
ship sail, windmill arm, etc. We are interested here in maximizing the pres- 
sure force of the wind consisting of non-interacting particles. An example of 
such a wind is provided by the flow of solar photons incident on a solar sail. 
We first studied the problem numerically and in a restricted class of bodies, 
and then, basing on the gained experience and intuition, found the solution 
analytically. 

On the picture below, there are shown "pre-optimal" bodies, that is, 
elements of maximizing sequences of bodies for problems (A) and (B). (In 
the case (B), K is the unit circle.) 
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(a) (b) 

One can see that the boundary of the second figure is much more compli- 
cated than of the first one. "Complexity" of the boundary can be measured 
by the maximal value of rotation angle of the normal vector on a small part 
of the boundary. This "maximal instantaneous rotation" is, approximately, 
90° + ~ 129° for the first figure and 360° for the second one. 

To be precise, fix a convex bounded body K and consider a sequence 
B n , n = 1, 2, . . . approximating K. Define the value rot({-B n }) measuring 
boundary complexity of the sequence in the following way. For each point 
x G dB n , let v™ be the outer unit normal at x and let 9 n {x) := axg(v™\. 
For any x G dK, define rot n (x,e) := sup xl>X2edBnfWsix) (6 n {x 2 ) - n (zi))B 
let us call it (x, e)-rotation of the boundary dB n . It is monotone non- 
increasing as £ —* 0. Then put rot n (e) = sup^^ rot n (x, e) and define 
rot(e) = limsupjj^oo rot n (e): e-rotation for the approximating sequence 
B n . Finally, define rot({-B n }) = lim £ ^o r ot(e) and call this value rotation 
of boundary for the given sequence of bodies. This value is tc/2 + \l/ f° r the 
first sequence of figures, and 2n for the second one. The difference is almost 
threefold. 

Further, the boundary length for the second sequence of figures tends to 
infinity. On the other hand, one easily calculates that the boundary length 
tends to cos ^/ 2 ) = 1-501 • 2ir for the first sequence; that is, the limit value 
is approximately 1.5 times the perimeter of unit circle. 

2 Recall that Af £ (x) is the e- neighborhood of x. 
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Throughout this paper it was assumed that the "inner temperature" of 
the wind is zero. However, one can show that in the case of positive temper- 
ature, that is, chaotic relative motion of wind particles, the functional to be 
minimized is proportional to the functional JF (J5J). Therefore all the results 
of this paper remain valid in the case of positive temperature; the detailed 
explication is postponed to a future paper. We are also going to study the 
three-dimensional case in a future paper. 



Appendix 1 

Denote by B^ the set depicted on Fig. (a) in Introduction. Here e divides n 
and conv-B^ is a (7r/e)-sided regular polygon inscribed in the unit circle K\. 
Its perimeter is ^(convB^)] = 2n sine/e. The convex part of dB^ is a fi- 
nite collection of points and all the cavities of B^ e > are isosceles right triangles 
similar to Q — {(xi, x 2 ) : < Xx < 1, < x 2 < 1/2 — |l/2 — Xi\}. Applying 
®, one gets K(B&) = -\d(oonvB&)\ ■ ■ e 2 = -2tt (sins/e) .F(fi) ■ e 2 . 

According to the same formula ([6]), resistance of the unit circle K\ equals 
IZ(Ki) = —\dKi\ ■ e 2 = — 2tc • e 2 . It remains to calculate jF(fi). 

Let A±, A 2 , and A\ 2 be the subsets of [—tt/2, tt/2] x [0, 1] given by the 
inequalities £ < — tamp, £ > 1 — tamp, and — tamp < £ < 1 — tamp, 
respectively. One easily verifies the following. 

(a) If (tp, £) G A\ then the corresponding billiard particle makes only one 
reflection from the left leg of Q and (pa(ip,^) = —tt/2 — (p. 

(b) If G A 2 then there is a single reflection from the right leg of Q 
and <pn(<p,g) = tt/2 - (p. 

(c) If (<£>,£) G Ayi then there is a double reflection, either from the left 
and then from the right leg or vice versa, and <pn(<p, £) = (p. 

Therefore, according to (151) , one has 

jr(fi) = - if (1 + C os(vr/2 + 2^)) cosy? dpd£+ 
+ - // (1 + cos(2<y? — 7r/2)) cos (p dip d£ + — if 2 cos ip dip dC, = 

8 J JA 2 8 </ ^Aia 

= I + 11 + III. (8) 
Direct calculation gives I — II — | — /// = | (\/2 — 1); thus, jF(fi) = 
V2. 



18 



Appendix 2 



Recall that a right triangle is called canonical if (a) it is situated above its 
hypotenuse and (b) the median dropped on the hypotenuse is vertical. The 
angle a the hypotenuse forms with the horizontal line, a G (— n/2, 7r/2), 
is called inclination of the triangle. Consider a particle that intersects the 
hypotenuse, gets into the triangle, makes one or two reflections from the legs, 
and then intersects the hypotenuse again and leaves the triangle. Denote by 
ip the angle the initial velocity v forms with the vector e 2 = (0, 1), and by 
ip + , the angle the final velocity v + forms with — e 2 . Thus, one has v = 
(sin <f, cos (p) and v + = — (sin p + , cos ip + ), where f and ip + vary between 
— 7r/2 + a and ir/2 + a. 

Parametrize the hypotenuse by the parameter £ G [0, 1]; the value £ = 
corresponds to the left endpoint of the hypotenuse, and the value £ = 
1, to the right one. Like in Appendix 1, denote by A x the set of values 
(</?, G [— tt/2 + a, tt/2 + a] x [0, 1] corresponding to particles having a 
single reflections from the left leg, by A 2 , the set of values corresponding to 
a single reflection from the right leg, and by A\ 2) the set corresponding to 
particles having double reflections. One easily finds that A\ is given by the 
inequality £ < A 2 , by the inequality £ > 1 - ^g^y, and A 12 , by 

the double inequality - co ^ a) < £ < 1 - co ^ a) ■ Moreover, for (<p,£) G A 1 
holds ip + = — 7r/2 + a — (p, for (</?,£) G A 2 holds y? + = 7r/2 + a — <p, and 
for (</?, G Ai 2 holds v? + = (p. On the picture below, ip < 0, ct > and 




Consider the parallel beam of particles falling on the hypotenuse in the 
direction (p. If | sin</?| < cos(</? — a), then the portion of particles that make 
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only one reflection equals | sin</?|/ cos(y? — a) and the direction of reflected 
particles is ±7r/2 + a — ip; one has to choose the sign "+" if sirup > 0, and 
"— " if sirup < 0. The rest of particles make double reflections; the portion 
of these particles is 1 — |sin</?|/ cos(y? — a), and the direction of reflected 
particles is <p. If | sin 1 > cos(p — a) then all the particles make a single 
reflection and the direction of reflected particles is ±n/2 + a — p. 

Now, consider the arc of circumference of angular size 2\1/ contained in the 
half-plane X2 > 0, with the endpoints A = (0, 0) and C = (1, 0). Parametrize 
this arc with the parameter a e [— the value a = — ^ corresponds to 
the point A, and a — to the point C . Divide it into a large number of 
small arcs and substitute each of them with two legs of the corresponding 
canonical triangle. The resulting broken line (shown on the figure below) 
defines a standard cavity. Denote by 5 the maximum length of a small arc. 



B 




For small 5, the scheme of billiard reflection can be approximately sub- 
stituted with the following description (pseudo-billiard reflections from the 
arc ABC). A particle of some mass moving in a direction p e (— 7r/2, 7r/2) 
is reflected from the arc ABC. If | sin</?| < cos(p — a), it is split into two 
"splinters" of relative masses | sin <^|/ cos(p — a) and 1 — | sin</?|/ cos(p — a). 
The first splinter is reflected in the direction ±7r/2 + a — <p, and the second, 
in the direction tp. If | sin <y?| > cos(p — a), there is no splitting, and the 
whole particle is reflected in the direction ±7r/2 + a — p. The described 
dynamics will be called pseudo-billiard one. A particle of unit mass starts 
moving at a point of X = AC in a direction ip e (— n/2, vr/2), and after 
several pseudo-billiard reflections, the resulting splinters return to X. 

As a result of the described substitution of the billiard dynamics with the 
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pseudo-billiard one, one obtains the function 

° i J-tt/2J0 

where = rrii((p,C) are masses and <pf = (pf(<p,£), final directions of the 
splinters resulting from the particle with the initial data (<£>,£). (As we 
will see later, splitting can really occur only once, after the first reflection, 
therefore there are at most two splinters.) The difference between Ffif) and 
the true value of the functional T Q5J) is 0(S), 5 — ► + . This fact can be 
expressed as lim,5^o+ ^(broken line) = J 7 ^). Below we will calculate 

In order to describe the pseudo-billiard motion, it is helpful to change the 
variables. Consider the circumference containing the arc under consideration, 
and parametrize it with the same angular variable a; this time a varies 
in [— 7r, 7t]. Consider a particle that starts moving at some point (3 of the 
circumference, intersects I at some point (£, 0), and then reflects from the arc, 
according to the pseudo-billiard rule, at a point a. Thus, one has \I> < \j3\ < n 
and |a| < ^. If \f3\ < n/2, there is no splitting, and if \j3\ > ir/2, there is. 

Let us describe the dynamics of the first splinter. For a while, change 
the notation; let (3 =: a_i, a =: a , and designate by ax the point of 
intersection of the splinter trajectory with the circumference. Denote by if 
the initial direction of the particle, and by ip', the direction of the splinter 
after the first reflection. (We do not call it ip + , since there may be more 
reflections.) One has if = ±n/2 + (a + a_i)/2 and ip' = tt/2 + (a + «i)/2. 
Then, taking into account that ip' = ±n/2 + a — ip, one gets 

ip = (a — «i)/2 and (a_i + «i)/2 = 7r/2, 

the equalities being true mod7r. In other words, the points a_i and a.\ lie 
on the same vertical line; see the figure. 
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If cci belongs to the arc [— ty, ^] then there occurs one more reflec- 
tion, this time without splitting, since the splinter arrived from the point 
a G [— 7r/2, 7r/2]. Extend the trajectory after the second reflection until 
the intersection with the circumference at a point «2- Using an argument 
analogous to the one stated above, one derives the formula ao + «2 = n', it 
follows that the point «2 does not belong to the arc, that is, there are no 
reflections anymore. 

Summarizing, the pseudo-billiard dynamics is as follows. After the first 
reflection from the arc, the particle may, and may not, split into two "splin- 
ters". If a_i G [— 7r/2, — \&] U [ty, 7r/2], there are no splitting, and the 
reflection is unique. If a_i G [— it, —ir/2) U (n/2, n], there is splitting into 
two splinters. If a_x G [— n + ^, —n/2) U (7r/2, 7r — \&], the first splinter 
makes no reflections anymore. If a_i G [— 7r, — 7r + ^] U [7r — ^, 7r], it makes 
one more reflection (without splitting) from the arc, and the final direction 
is ip + = ti/2 + («i + «2)/2. Taking into account the above equalities, one 
gets ip — tp + = a:_i + ao + n. 

Note that the factor 1 + cos(v? — y? + ), meaning the impact force per unit 
mass, equals 2 for the second splinter. For the first splinter that makes no 
reflections, as well as for the reflection without splitting, this factor equals 
1 + cos((^ — ip') = 1 + | sina:_i|. Finally, for the first splinter that makes one 
more reflection, this factor equals 1 + cos(y? — tp + ) = 1 — cos(a;o + 

Let us pass from the variables ip and £ to a = a and f3 = a_i and 
calculate the integral J 7 ^) in terms of the new variables. The points a 
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and 8 on the circumference have the cartesian coordinates n 1 , T , (sin \I> + 
sin a, — cos \& + cos a) and 2 ^ (sin \1> + sin 8, — cos \I/ + cos /3) , respectively. 
The interval with the endpoints a and 3 intersects with the interval X at the 
point (£, 0), where 

^ _ sin(# + g) - sin(^ + 3) + sin(/3 - a) ^ 
2 sin \1/ (cos a — cos /3) 

Further, one has 

a + 8±u 

V = 2 ' ^ ^ 

one has to take the sign "— " or "+", if 3 > or 3 < 0, respectively. 
Therefore, cos<£> = | sin ^r-\- 

The point (a, (3) runs the set [-^, ^] x ([-vr, U it]), and the 
mapping (a,/3) ^ (y?,0 given by fl9l). f|T0l) is a one-to-one mapping from this 
set to [— 7r/2, 7r/2] x [0, 1], with the Jacobian 

D(<p, cos(^ + a) + cos(^ + /5) - 2 cos(/9 - a) 
D(a,j3) 4 sin \l/ (cos /3 — cos a) 

sin(^ + /5) - sin(^ + a) + sin(a - 8) . 

+ — , • t / a \T — sma + sm/3 = 

4 sm V (cos p — cos ay 



sin^ 



4sin^ sin^ ' 



(11) 



this implies that the integration factor equals 

3 — a 



cos (p dip 



sm 



da dd. 



4sin\I> 

Further, the mass of the first splinter is | cos \ j | sin | , and of the second 

one, 1 — | cos j | sin Note also that integrating over 3 G [— 7r, — >&] U 

[iff, 7r] can be substituted with integrating over 8 G [iff, ir] with subsequent 
duplication of the result. With this substitution, one always has sin > 0. 

The integral J 7 ^) can be written down as the sum J 7 ^) = I + II + 
77/ + IV, where 



16 sin 



— — / da I (1 + snip) sin dp, 
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II 



16 sin * 



01 • $ ~ a 

2 sin 



a + (3 



III 



IV 



16 sin* J_^ 



da 



da I 'l | sm — ; cos 

* Jtt/2 

>!' /• 71 ^ _|_ o 

da / (1 + sin/?) cos — - — d(3 



dp, 



tt/2 

(1 - cos(a + (3)) 



16 sin* J_m 
As a result of simple calculation, one obtains 



cos ■ 



a + f3 



dp. 



I = III 



16 sin * 



4 sin* 



1V2 . * 16 . 4 * 



16 sin * 



sin sin — 

3 2 3 2 



16 y/2 sin — - 8 * 
2 



IV 



- sin* + - sm d * + - * 

3 9 3 



16 sin * 

Summing these expressions, one finally comes to the formula (JTj): 

/ y \ , 1 • 2 t 2^ sin f - 2 sin 4 f - * 

jr(^) = 1 + - sin 2 * + ? 2 _ 

v 7 6 sin* 
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